Phase diagram of chiral quark matter: Fulde-Ferrell pairing from weak to strong 
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We calculate the phase diagram of two-flavor quark matter in the temperature-flavor asymme- 
try plane in the case where there are three competing phases: the homogeneous Bardeen-Cooper- 
Schrieffer (BCS) phase, the unpaired phase, and a phase with broken spatial symmetry, which is here 
taken to be the counterpart of the Fulde-Ferrell (FF) phase in condensed matter physics. The sys- 
tem belongs to the universality class of paramagnetic-ferromagnetic-helical systems, and therefore 
contains a tricritical Lifshitz point in its phase diagram, where the momentum scale characterizing 
the breaking of translational invariance has a critical exponent of 1/2 to leading order. We vary the 
coupling constant of the theory, which is obtained from integrating out the gluonic degrees of free- 
dom. In weak coupling, the FF phase is favored at arbitrary flavor asymmetries for sufficiently low 
temperatures; at intermediate coupling its occupancy domain is shifted towards larger asymmetries. 
Strong coupling features a new regime of an inhomogeneous FF state, which we identify with a 
current-carrying Bose-Einstein condensate of tightly bound up and down quarks. The temperature 
and asymmetry dependence of the gap function is studied. It is shown that the anomalous tem- 
perature dependence of the gap in the homogeneous, flavor-asymmetric phase is transformed into 
a normal dependence (self-similar to the BCS phase) at arbitrary coupling, once the FF phase is 
allowed for. We analyze the occupation numbers and the Cooper-pair wave function and show that 
when the condensate momentum is orthogonal to the particle momentum the minority component 
contains a blocking region (breach) around the Fermi sphere in the weak-coupling limit, which en- 
gulfs more low-momentum modes as the coupling is increased, and eventually leads to a topological 
change in strong coupling, where the minority Fermi sphere contains either two occupied strips or 
an empty sphere. For non-orthogonal momenta, the blocking region is either reduced or extinct, 
i.e., no topological changes are observed. 
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I. INTRODUCTION 

The phase diagram of dense hadronic matter is of great 
interest in relativistic astrophysics of compact stars and 
heavy- ion collisions. Because of the quark substructure of 
nucleons predicted by quantum chromodynamics (QCD), 
nuclear matter will undergo a phase transition to quark 
matter if squeezed to sufficiently high densities. In the 
quark matter phase the "liberated" quarks occupy con- 
tinuum states which, in the low-temperature and high- 
density regime, arrange themselves in a Fermi sphere. In 
this regime quark matter will exhibit the phenomenon of 
pairing, or color superconductivity, because the quark- 
quark interaction is attractive in certain channels (early 
work on this subject is reviewed in Ref. p|; various recent 
aspects are reviewed, e.g., in Refs. @|). 

Experimental verification of deconfined quark matter 
and color superconductivity can be found in the phe- 
nomenology of compact (neutron) stars 0, 0, H, Is 0, 
H, H EH QJj • The remnants of core-collapse supernovae 
is a hot and dense core of hadronic matter (the proto- 
neutron star), which rapidly cools from temperatures of 
the order of 100 MeV down to temperatures of the or- 
der of 1 MeV. The matter (in the baryonic or the de- 
confined quark matter phase) is initially nearly isospin 
symmetric. Electron capture on protons (it-quarks in de- 
confined quark matter) drives the isospin chemical poten- 



tial up to values of the order of 100 MeV. Consequently, 
in parallel to rapid cooling, the isospin asymmetry of 
matter grows rapidly. In effect, compact stars traverse 
the temperature-isospin chemical potential plane from 
the high-temperature, low-isospin domain down to the 
low-temperature, high-isospin domain. One of the prime 
interests of this work is to compute the phase diagram 
of superconducting quark matter within this plane. The 
central problem is, of course, to predict the asymptotic 
state of matter reached during this evolution, since cur- 
rently observable compact stars are in this asymptotic 
state, where they cool slowly from temperatures of the 
order of 100 keV down to 10 keV with nearly frozen com- 
position [II El El El Es 1 - 

The quark-quark interaction is strongest for pairing 
between up and down quarks which are antisymmetric 
in color [H, Q . At intermediate densities quark matter is 
composed of up and down quarks, while strange quarks 
can appear in substantial amounts at higher densities. 
As the isospin chemical potential grows and reaches the 
asymptotic value determined by the /?-equilibrium con- 
dition /Xd — — n e among d and u quarks and electrons, 
where fii with i = d, u, e are the chemical potentials, the 
Fermi spheres of up and down quarks are shifted apart. 
The difference can reach values ~ 100 MeV, the magni- 
tude of the electron chemical potential in /3-equilibrated 
matter. Thus, the cross-flavor pairing must overcome the 
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disruptive effect of mismatched Fermi surfaces. Further- 
more, if the physical strange quark mass is close to its 
current mass M s ~ 100 MeV, electrons may be gradually 
replaced by strange quarks, which again will disfavor the 
cross-flavor pairing. Such a situation can be described by 
assuming that the densities (or alternatively the chemi- 
cal potentials) of up and down quarks are unequal, as we 
shall do throughout this paper. 

A superconducting phase with asymmetric cross- 
species pairing needs to optimize the overlap between the 
Fermi surfaces (which is perfect in the symmetric BCS 
state) to attain the maximum possible condensation en- 
ergy. If temperature and correlations are negligible, an 
infinitesimal shift in the spherically symmetric, concen- 
tric Fermi surfaces is sufficient to destroy the supercon- 
ducting phase. A less restrictive (and less symmetric) 
phase space configuration in which the condensate car- 
ries non-zero momentum with respect to some fixed frame 
is more flexible; the penalty in the energy budget for the 
(always positive) kinetic energy of condensate motion is 
compensated by the gain in the (negative) condensation 
energy. This class of phases was first studied by Fulde 
and Ferrell (FF) [l7|, who suggested a condensate order 
parameter that varies as A ~ exp(iq • r), and by Larkin 
and Ovchinnikov (LO) [18| . who explored a number of 
spatially varying parameters. The simple form of the 
FF order parameter permits to study this phase at ar- 
bitrary temperatures [19( . Larkin and Ovchinnikov used 
instead the Ginzburg-Landau theory, valid in the vicinity 
of a second-order phase transition, to analyze more com- 
plicated patterns of symmetry breaking (lattices). The 
quark matter counterparts of the FF phase have been 
discussed in Refs. [H E0, H H [H H| . The quark mat- 
ter counterparts of the Larkin-Ovchinnikov phase, within 
the Ginzburg-Landau approximation, have been studied 
in the two- and three-flavor cases in Refs. [H, HE HE HE 
l30l . l3lj . In three-flavor quark matter the face-center cu- 
bic (fee) lattices were identified as having particularly low 
free energy [28| in the regime where the Ginzburg-Landau 
(GL) theory is applicable. The spatial form of the con- 
densate beyond the GL regime is not known. Further- 
more, if the dynamics of the gluon field is retained, there 
appears the possibility of gluon condensation instead of 
the onset of macroscopic supercurrents [H, [32], HE H3 . 
A more symmetric superconducting phase exploits defor- 
mations of Fermi spheres as the mechanism of restoring 
the coherence needed for pairing at the cost of kinetic 
energy loss caused by these deformations. The deformed 
Fermi surface phase requires minimal breaking of spa- 
tial symmetry and is more symmetric than the lattice 
phases above [H, HE HE HI] . To lowest order the defor- 
mations of the Fermi spheres break the rotational 0(3) 
symmetry down to 0(2). In this paper we shall adopt a 
particularly simple form of the order parameter: the sin- 
gle plane-wave FF ansatz. It shares many features with 
the more complex phases mentioned above and permits 
us to study the phases with broken spatial symmetry in 
a straightforward manner. 



Our second key interest in this paper is to explore 
the phase diagram as the coupling is increased from 
its conventional value to larger values. To our knowl- 
edge this is the first investigation of asymmteric, in- 
homogeneous quark matter at strong coupling. As is 
well-known, for non-relativistic fermions the Nozieres- 
Schmitt-Rink [3^| theory describes the mean-field evo- 
lution of the phase diagram from the BCS phase to the 
Bose-Einstein condensate (BEC) regime. This theory has 
been generalized to the case of unequal populations of 
fermions [40| . Furthermore, there are a number relativis- 
tic treatments of the problem for both matching and mis- 
matched Fermi surfaces [HI, H EE H EE EE EE EH • 
Here we wish to extend the emerging picture of of the 
phase diagram by including inhomogeneous, relativis- 
tic, phases of superconducting chiral fermions (quarks). 
While weak-coupling studies are likely to be sufficient 
for neutron star phenomenology, other superconducting 
systems may help to explore (experimentally) the strong- 
coupling regime as well. Examples are pairing in dilute 
alkali gases [M HE HE HE HE 0, HE HE HE HI , where 
density-imbalanced systems were created, for example, 
by using two hyperfine states of 6 Li fermionic isotopes, 
and graphene where the quasiparticle spectrum has rela- 
tivistic form [HE HE HE H^] • 

Color superconductivity of quark matter can find its 
verification in the astrophysics of compact (neutron) 
stars and supernovae. The phase with broken spatial 
symmetry, such as the FF phase studied here or the 
crystalline phases with complex lattice structure, have 
their specific signatures which we shortly describe below. 
One avenue of discerning these phases in compact stars 
is through their cooling behavior [IE HE HE HE HE HI] • 
Normal quark matter cools too fast via the direct URCA 
process to account for measured neutron star surface 
temperatures. This implies that either quark stars are 
unobservable in the X-ray spectrum or quark matter is 
paired, the pairing gap preventing rapid cooling. Since 
the gaps are much larger than the relevant temperatures, 
pairing will effectively block the relevant neutrino pro- 
cesses. However, for the phases with broken spatial sym- 
metry, not the entire Fermi surface is "gapped" ; if this is 
so, cooling will take place at a reduced rate from the seg- 
ments of the Fermi surface that are ungapped. In fact, 
it is possible to interpolate between the fully gapped and 
unpaired limits by assigning some amount of "gapless- 
ness" in the spectrum of quasiparticles [6J, l65l |66( . An- 
other avenue is offered by the fact that crystalline color 
superconducting (CCS) matter may have shear moduli 
that are by orders of magnitude larger than those in 
neutron star crusts [69| . Consequently, CCS matter can 
support quadrupole deformations and can be a source of 
gravitational waves [jE HE [72[ . Upper limits on gravi- 
tational waves derived for known pulsars are within the 
range where the hypothesis of CCS matter in the cores 
of neutron stars can be tested [73| . 

Ab initio methods, such as lattice QCD, are (cur- 
rently) not useful for studying the quark matter phase 
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diagram at high densities. Furthermore, in the den- 
sity regime between twice and ten-fold nuclear satura- 
tion density, the coupling is too strong for perturbative 
QCD to be quantitatively accurate. Qualitative predic- 
tions are based either on the single-gluon exchange model 
or point-coupling effective models, such as the Nambu- 
Jona-Lasinio (NJL) model Q. In this paper we shall 
start with the QCD Lagrangian, which contains explicit 
gluon fields and derive the Dyson-Schwinger equations 
describing the quark and gluon fields, and only at a later 
stage we will integrate them out. Such a methodology has 
been applied to homogeneous superconductors in Ref. Q . 
Clearly, our final equations could have been also obtained 
directly from the NJL model Lagrangian without explicit 
reference to the QCD Lagrangian. 

This paper is organized as follows. In Section [II] we set 
up the formalism to treat the FF and related phases in 
chiral quark matter. We start with the QCD partition 
function, perform a local transformation on the fields to 
allow for motion of pairs with non-zero momentum. The 
pairing fields are introduced via bosonization of the ac- 
tion after which we take the mean-field approximation. 
We write down the general Dyson-Schwinger equations 
for the quark and gluon propagators for our systems, af- 
ter which the gluon propagator is postulated to have a 
contact form, which in turn eliminates the gluon fields 
and reduces the model to the four-fermion coupling NJL 
model. Our final equations (f42)) - (|48f constitute four in- 
tegral equation that must be solved simultaneously. Sec- 
tion [TTT] discusses the numerical results for the phase di- 
agram, gap function, occupation numbers, and the ker- 
nel of the gap equation with attention to the differences 
between the strong- and weak-coupling regimes. In par- 
ticular, we study the phase-space topological structure 
of the minority component in three interesting cases of 
weak, intermediate and strong coupling. We close with a 
discussion of our results in Section IIVI 

Appendix lAl shows the details of the transformation of 
the characteristic equation for the quasiparticle modes in 
the FF superconductor from the form that contain the 
propagators of the theory to a form that contains disper- 
sion relations in the normal state and the gap function. 



II. FORMALISM 

The theory is described by the following functional in- 
tegral 



where 



Cqcd = ^(iVfy - m)V - -(Fp 2 + g^A^, (2) 

and i\) is the quark field, = A^X a /2 is the gluon field, 
A a are the Gell-Mann matrices, g is the strong coupling 
constant, m is the quark mass, F" = d^A^ — d v A^ + 
g pbc A b^ A c ig the field strength tensor of the Yang-Mills 

field, f abc are the structure constants. To eliminate glu- 
ons we write the generating functional {1} 



Zqcd = I VipVipexp 
where 



d^xi/j^d^ -m)i/j + r\j] 



(3) 



r[j] =log J VAexpfj d 



' (4) 

and j% = ■0T a 7 M -0 is the color current of quarks, T a = 
A a /2 are the SU(3) generators. Expanding the action ((H) 
in powers of the quark current and keeping the leading 
order quadratic term we obtain [74j 

r\j] = ^J d^xefiy jZ(x)D%(x,y)j h v (y), (5) 

where D?^(x,y) is the single-particle irreducible gluon 
correlation function in the pure Yang-Mills theory 

Df v {x,y) = (A^x)A b u (y)) - (A{x)l)(A b M) (6) 

and the brackets stand for average over the gluon field, 
i.e., 



JVAe^...(-lJF 2 ) 
JVAexp(-{jF*) ' 



(7) 



The partition function of the theory where the gluon field 
is integrated out is given by a path integral over the quark 
fields tp(x) as 



Z= P^XVexp {S $, V>] + Si[i>^]} 



(8) 



AQCD — / V^V^VA a exp / d^xCqcD , (1) 



where the first term in the action is the free Fermi-gas 
contribution and the second term is the four-fermion 
gluon-mediated interaction contribution 



S $,il>] = I d A x d 4 yj)(x)[G+] 1 (x,y)i>(y), 



(9) 
(10) 
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where G^ix^y) is the free quark propagator, = T a j^ is the quark-gluon vertex, and ± superscripts on the quark 
propagator refer to fermions and charge-conjugate fermions. We consider hereafter the case of Nf — 2 flavors (u and 
d quarks) with iV c = 3 colors. Our basis is defined as 



1> 



(11) 



Consider a local transformation on the quark fields given by 



tp' —>■ ipe 



i8(x) 



ip — ► ipe 



\6(x) 



(12) 



We shall focus on the special case 6(x) = Q )Jl x fJ, /2, where the four vector Q = (0, Q) has only a non-vanishing spatial 
component. The transformation (|12p breaks the translational symmetry of the theory. The particular form of the 
transformation leads to the free quark propagator with non- vanishing momentum Q/2 with respect to some fixed 
frame. Under the transformation (|12p the kinetic part of the partition function retains its form, however the free 
quark Green's function becomes (flavor and color indices are suppressed) 



(x, y) = [G+]- l {x,y) + [7^%)] S(x - y). 



(13) 



The effect of the transformation becomes transparent if we examine the Fourier transform of the propagator (|13p 



\G 



±1-1 



7^ (±Q M /2 + kft) ± M7o - m, 



(14) 



which is seen to be identical to the Fourier transform of the propagator [G^]" 1 with respect to the relative and 
center-of-mass space-time coordinates in a frame with non- vanishing center-of-mass momentum Q. Here /x is the 
chemical potential associated with the conserved baryon charge. The action is bosonized by introducing pair fields 



E + (x,y) = g 2 T^ c (x)$(y))nD a J> u (x,y), 



s (y,x) = 7o[A+(a;, y)] t 7o, 



(15) 



where f £ = G(r^) T G _1 = —T^j^. The interaction term of the action will be evaluated in the mean- field approxi- 
mation by keeping the leading-order term in the expansion with respect to the deviations of the two-point function 
K(x,y) = ipc(x)ip(y) from its mean value (K(x,y)}. In this approximation the partition function is given by 



Zmf = jvipVipexp{SE.[ip,i;]} 



x exp< 
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dxdyTT[ 7Q (K^y,x)) 7 o^(K(x iy ))^ b + {K(x,y))T^ (K^y,x)) l0 r^D a 'l(x,y) , (16) 



where the trace is over Dirac, color, and flavor indices and 



SsbP,ip]= / dxdyl ip{x)[G£] 1 {x, y)ip{y) + -[ip c (x)E+(x, y)ip(y) + ip(y)Z (y,x)ip c {x)] 



We now rewrite the functional integral in terms of Nambu-Gorkov spinor fields 



* = 



ipe 



* = (ip,rp c ) 



(17) 



(18) 



where %pc = Cip T is the charge-conjugate spinor; C is the charge-conjugation matrix. In this basis, the full propagator 
and the self-energy are given by 



G+ F~ 
F+ G' 



5+ S- 



(19) 
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Their elements are related by the Schwinger-Dyson equations 



\G d 



Gn 



F ± = - 



G 



Gj 



(20) 
(21) 



where is the so-called anomalous propagator. We will ignore the normal self-energy in the following. At the 
mean-field level the partition function is then given by 

£ MF = [det.^a-^^expj^ J j -^Tr [F~ (k)f^F+ (p)rfl D%{k - p)} , (22) 

where /3 is the inverse temperature, d 4 /c = dkod 3 k, where the ko integration is understood as summation over the 
discrete fermionic Matsubara frequencies; note that the trace runs over color, Dirac, flavor spaces; the trace over the 
Nambu-Gorkov space has been carried out. 

In the following we shall restrict our discussion to the chiral limit of massless quarks. In this case the free quark 
propagator is diagonal in color and flavor space 



G = dia; 



(/^f±ii ✓"fieZ >^f±u f-t±d fi±u /^i±d\ 



(23) 



The gap functions, in the basis (HJ), are given by 
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V 














0/ 



(24) 



We shall expand the gap matrices in the basis of positive and negative states 

S+ 2 (fc)=]TAt 2 (fc)A e (fc), 



and negative e 



(25) 

energy states in 



where A e (fc) = (1 + ejoj ■ k/\k\)/2 are the projectors onto the positive e = 
the chiral limit. 

The full G propagator is diagonal in color-flavor space and its elements can be read off from the Schwinger-Dyson 
Eq. (|20|) . The full anomalous propagator has a matrix form identical to Eq. ([24| 



F ± = 



with the elements that follow from the Schwinger-Dyson equation ([21 



/ 











rp±ud 
rg 


0\ 










rp±du 
rg 















rp±ud 

gr 










j^±du 

gr 





























V 














0/ 



(26) 



±ud 
rg 

±ud 
gr 



/-f^u^i^Vifi rp±du _/-i z fd'^±s : y±u 

~~ Lt 0o ^2 •> r qr ~ U 0q ^1 L7 r ■ 



yt±du /~i^fd^±s : y±U 

rg — ^Or 1 - J 2 .(? ' 



(27) 
(28) 



The quasiparticle excitation spectrum is given by the eigenvalues of the fermionic determinant in Eq. (|22p . Equiva- 
lently, it is given by the poles of the propagators in Eqs. (|2"0"|) and (|2"Tj) , i.e., by [Cr] -1 = [f 1± ] _1 = 0. For the paired 
quarks this translates into the characteristic equation for the modes 



r<+9 



-£|A e ^| 2 A e (fc) = 0, 



(29) 
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where the indices i,j and /, g label the (admissible combinations of) color and flavor indices, respectively. The solution 
to Eq. (|2"5|) reads (see Appendix [A] for details) 



(30) 



where Es, e and Ea, £ are the parts of the spectrum which are even (symmetric) and odd (asymmetric) under exchange 
of quark flavors in a Cooper pair. These are given by 

,2 



£s, e (|fc|,|Q|,M) 2 = (\k\-ep,Y 
EA,e(\Q\,0,5n) = 5v + e\Q\cos9, 



(31) 
(32) 



where 9 is the angle formed by vectors k and Q, Sfi = (/i,; — (J>j)/2 and p, = (fa + fij)/2. The spectrum of unpaired blue 
quarks can be obtained as above by neglecting the contribution from the anomalous self-energy to the propagator. 
We quote only the final result, which is given by Mb) — Ef e (k,0,0, fit), where fib is the chemical potential of 

blue quarks. 

In the following we reduce our model to the NJL model featuring a contact four-fermion interaction. The momentum- 
space gluon propagator for a contact interaction has the simple form 



7->ah _ cab 9jiU 



(33) 



where A is a characteristic momentum scale. Upon inserting this result into Eq. (|22|) . evaluating the fcrmionic 
determinant, and carrying out the summation over the Matsubara frequencies we find 



IiiZmf = hiZ: 



A 

MF 



(34) 



where the first term includes the contribution from the condensate of red-green quarks and the second term is the 
contribution from the non-interacting blue quarks. The first term is given by 

lnZ^ F = ^£(|A+| 2 + 3 £ A+A+) 

+ \T,J ^ {P [E+ (K) ~ E~ (A^)] - 2 In f [-E+ (A* n )]- 2 In f[-E^ (A*)]}, (35) 

where f(uj) = [1 + exp(/3w 

\nZ° MF = 2 J ^L]T{_ln/ [-£+(*:, Mb,/)] ~ W K'^W,/)]} 



is the Fermi distribution. The contribution of the blue quarks is 
d 3 k 



(36) 



The thermodynamic potential is obtained from the logarithm of the partition function as 

1 



MF 







-In Z 



MF- 



(37) 



The stationary point (s) of the thermodynamic potential determine the equilibrium values of the order parameters 



MF 



= o, 



MF 



= 0, 



MF 



= 0. 



(38) 



dA\ dA e 2 a\Q\ 

The direction of the vector Q is chosen by the superconductor spontaneously. We shall work at fixed baryonic density 
and nonzero temperature. We shall consider variations in the relative concentrations of the u and d quarks, which are 
determined from the derivative of the thermodynamic potential with respect to the corresponding chemical potential 
as 

<9J7mf dflMF 



The explicit form of the gap equations follows upon minimization as 



Ai + 3A 2 = - 
3A X + A 2 = - 



2.9 



E 

e 

; U 2 Zl 



3A 2 
2.9 2 



d 3 p 



2A| 



(27r)3s+(Af)-S e -(Af) 



d 3 p 



2A% 



(2tt) 3 Et{Al)-Ee{Al) 



tanh 



tanh 



EtiAl) 



tanh 



tanh 



E-(At) 
E-(Al) 



(39) 

(40) 
(41) 



7 



Note that the minima are degenerate with respect to the particle charge index e, consequently we have dropped this 
index on the left-hand sides of the gap equations. Because the kernels on the right-hand side of the two gap equations 
are invariant under the interchange 1 ^> 2 the two gap functions must be degenerate Ai = — A2 = A, i.e. we are left 
with a single gap equation 



A 



3A 2 



E 



d 3 P 



(2nrE+(A-)-Ee(A-) 



tanh 







£+(A e ) 



tanh 



E~{A e ) 



(42) 



In the symmetric limit _E+(A e ) = E e (A e ) and when \Q\ = the gap equation reduces to the ordinary BCS gap 
equation for ultrarelativistic fermions. The densities of paired up and down quarks are given by 



d 3 k fdE+(A° n ) 



Ud 2 2^ J (2tt)3\ dn d 

e,n v 



tanh 







n b 1 ^ 



2^ J (^F\^- tanh 



where rib is the number density of blue quarks and 

dEt(Al 2 ) _ 

and 



it 



eE s , 



tanh 



tanh 



^E-(Al) 

~Ee(A* n ) 



Eh 



|A li2 | 2 




dEt(Al 2 ) 

The contribution of the unpaired quarks to the density is given by 

d 3 k 



n b = 



2 J 7^3E{/Ke + (^M b ,/)]-/[C(fc,^/]}- 



(43) 
(44) 



(45) 



(46) 



(47) 



The total quark density is n u + rid + n b = 3ns, where tib is the baryonic density. Finally, the equation for total 
momentum is given by 



E 



d 3 k 



ee cost 



Es,< 



dE s , e 



Eh 



d\Q\ 



tanh 



E% (A e ) 



(48) 



The zero temperature counterparts of the Eqs. (|42)) - (148|) are straightforward to obtain with the help of the identity 
tanh(cc/2) = 1 — 2/(x) and the limiting from of the Fermi distribution function f(x) = 0(—x) as [3 — > 00, where 9(x) 
is the Heaviside step function. 

Massive strange quarks can be present in matter, but may not participate in the pairing with the u and d quarks. 
This requires the strange quark mass to be larger than the chemical potential mismatch, since otherwise the cross 
flavor pairing with strange quarks can not be ignored (in other words, we ignore the possibility of the color- flavor- 
locking type pairing) . Under such circumstances the contribution of strange quarks to the thermodynamical potential 
is straightforward to write down 



(49) 



where E^ e = ^k 2 + m 2 . Here, m s is the strange quark mass, which is expected to lie within the range 100 < m s < 500 
MeV. ' ! 



III. RESULTS the four integral equations for the gap, the densities of up 



The model defined by Eqs. (|42|) - (|48f was solved numer- 
ically by finding a solution that satisfies simultaneously 
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and down quarks, and the nonzero momentum of the con- 
densate. The ultraviolet divergence of the integrals was 
regularized by a three-dimensional cut-off in momentum 
space \p\ < Ajvjl- The phenomenological value of the 
coupling constant G — g 2 /12A 2 in the (qq) Cooper chan- 
nel is related to the coupling constant in the (qq) di-quark 
Gd channel by the relation G = N C /(2N C — 2)Gd, where 
N c = 3 is the number of quark colors. The coupling 
constant Gd and the cut-off Ajvjx are fixed by adjusting 
the model to the vacuum mass and decay constant of the 
pion. We employ the parameter set Gd = 3. If GeV -2 
and G d A 2 NJL = 1.31. 

Our discussion below will be carried out at fixed bary- 
onic density ns = 0.55 fm -3 . We wish to explore the 
behavior of our model under variations of temperature, 
density asymmetry (or, equivalently, the asymmetry in 
the chemical potentials) and the magnitude of the cou- 
pling. Our numerical runs search for solutions in the 
density-asymmetry and temperature plane at some fixed 
value of the coupling constant. The results are shown for 
a coupling constant parameterized by the (dimensionless) 
ratio / of the actual coupling parameter to that in the 
weak-coupling theory, i.e., G. We consider three values of 
the dimensionless coupling / = 1.0, 2.0, and 2.93 which 
cover the range from weakly coupled theory (/ = 1) to 
the strong coupled theory / = 2.93. The latter value is 
slightly below the zero-temperature critical value of the 
coupling constant at which the FF phase is extinct. 



A. The phase diagram 

The phase diagram of the model in the density- 
asymmetry - temperature plane is shown in Fig. [1] for 
the (dimensionless) couplings / = 1.0, 2.0, and 2.93. The 
lowest value / = 1 corresponds to the weak-coupling the- 
ory where the occupation numbers of quarks have well de- 
fined fermionic nature (occupied states below the Fermi 
energy and empty states above). For the largest value 
/ = 2.93 these features vanish and the occupation num- 
bers are smooth functions of momentum. For obvious 
reasons, this limit will be referred to as the strong cou- 
pling limit. The phase diagram contains three phases: 
(1) the unpaired state, where the gap A and the mo- 
mentum Q vanish, (2) the BCS state, where A ^ 0, but 
Q = and (3) the FF state, where A/0 and Q ^ 0. 
These phases meet at a tricritical Lifshitz point (here- 
after L point) which divides the critical line of phase 
transitions to the normal state into two segments. On 
the small-a segment the normal state transforms to the 
BCS state, which is uniform in space (Q = 0). On the 
large-a segment the unpaired phase transforms into the 
superconducting phase with non-uniform order parame- 
ter (Q ^ 0). The pair- momentum goes to zero as the L 
point is approached along the critical line. An expansion 
of the free energy of the system for small order parameter 
near the L point contains terms F ~ ci|V^| 2 -|- C2 |v 2 -0| 4 . 
Since on one side of the L point the minimum of the free 




FIG. 1: (Color online) The phase diagram of chiral quark 
matter for three (dimensionless) couplings / = 1.0, 2.0 and 
2.93. For each coupling the x and y axis are normalized to 
the values of the critical asymmetry a c (f) at T = and the 
critical temperature T c (f), respectively. The values of these 
parameters are: a c (l) = 0.18, a a {2) = 0.77, a c (2.93) = 1.0 
and T c (l) = 23.0, T c (2) = 144.0 and T c (2.93) = 252.0 MeV. 
The upper left domains marked as BCS correspond to the 
BCS state with Q — 0. The lower right domains marked as FF 
correspond to the Fulde-Ferrell state with Q / 0. The upper 
right domain corresponds to the unpaired state ( A = 0) . The 
three phases meet at the (Lifshitz) tricritical point. 



energy corresponds to Q = 0, while on the other side 
the minimum corresponds to Q ^ the coefficient C2 
must change sign at the L point and C4 must be positive. 
This suggests that chiral quark matter belongs to the 
universality class of paramagnetic-ferromagnetic-helical 
systems for which these terms lead to the existence of 
the L point. The critical behavior of the vector Q was 
determined from renormalization group arguments near 
the L point [76| . Along the large-a branch the momen- 
tum of the condensate varies to leading order as Q ~ \a\@, 
where [3 = 1/2, a = (a — a c )/a c and a c is the critical 
value of the asymmetry. 

The structure of the phase diagram can be under- 
stood qualitatively in terms of the quasiparticle spec- 
trum, Eqs. (f3"Tj) and (f3"2")) . At nonzero a the Fermi sur- 
faces of the species are mismatched by an amount given 
by Eq. (f3"2"|) . In the homogeneous BCS state the mis- 
match is a scalar (5/i, independent of the direction in mo- 
mentum space. Because cos 9 assumes both positive and 
negative values, the mismatch is "modulated" in momen- 
tum space for the FF phase. One can speak of coherence 
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(or decoherence) to quantify to which extent the Fermi 
surfaces of the up and down quarks overlap. The coher- 
ence is maximal for perfectly overlapping Fermi surfaces 
(Ea — 0). While in the homogeneous phase the coher- 
ence decreases as 8/j, is increased, in the FF phase the 
decoherence can be compensated as the Fermi surfaces 
of the majority and minority components come closer to- 
gether in some directions, which in turn can increase the 
pairing energy. This, however, increases also the overall 
kinetic energy of the condensate. The compromise among 
these two effects can produce an overall gain in the neg- 
ative condensate energy, which eventually makes the FF 
phase more favorable than the BCS phase. Nonzero tem- 
perature "smears" the Fermi surfaces, i.e., it reduces the 
mismatch over the range \J3Ef(u, d)] , where Ef(u,cI) 
are the Fermi energies of the up and down quarks. 

Let us return the phase diagram and comment on its 
structure. The weakly coupled T — ► limit is charac- 
terized by the fact that the FF phase is favored at any 
asymmetry. At T = the excitations are confined near 
the Fermi surface and a small mismatch is sufficient to de- 
stroy the BCS state. Nonzero temperatures increase the 
Fermi surface smearing and the coherence, therefore the 
domain occupied by the BCS state is larger. Increasing 
the coupling has two effects: first, the BCS state becomes 
more robust, as the ratio of the gap A(T = 0) to the aver- 
age chemical potential p, increases; the occupancy domain 
of the FF state is correspondingly reduced. Secondly, 
the nature of the BCS state changes in the sense that 
the condensate evolves to a state that has characteristics 
of a homogeneous Bose-Einstein condensate (BEC). The 
nature of the FF state also changes; in the FF state ad- 
jacent to the homogeneous BEC state only the majority 
distribution has fermionic nature (a pronounced Fermi 
sphere); the minority component has a shell-structure 
distribution (for more details see the discussion below). 
It is evident that there exist a limiting coupling beyond 
which the FF state is vanishing. We note that a topo- 
logically similar phase diagram exists for non-relativistic 
fermions in weak coupling [77], [78[ • Ref. [z3] suggests an 
approximately self-similar evolution of the phase diagram 
in the strong-coupling regime, which is at variance to the 
results shown in Fig. [TJ 

If the transition to the color-superconducting state oc- 
curs at sufficiently high densities, newborn compact star 
will cross the diagram in Fig. [JJ starting from the upper 
left corner and move down the temperatures axis towards 
nonzero asymmetries. The key implication of the phase 
diagram is that the asymptotic low-temperature state 
corresponds to the FF state independent of the magni- 
tude of the flavor asymmetry. 

An alternative view on the phase diagram of paired 
chiral quark matter offers Fig. [21 where the phases are 
plotted in the temperature and chemical-potential mis- 
match plane. Consistent with what we have learned from 
Fig. [1] the domain of occupancy of the FF state shrinks 
as the coupling is increased, although there is no simple 
mapping between the two cases. For / = 1 and / = 2 
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FIG. 2: (Color online) The phase diagram of chiral quark mat- 
ter in the temperature - chemical potential mismatch (8fi/2) 
plane, for three couplings, / = 1 (upper panel), / = 2 (mid- 
dle panel), / = 2.93 (lower panel). In each panel the lower 
left domain corresponds to the BCS state with Q = 0, the 
domain between the solid and the dashed lines, corresponds 
to the Fulde-Ferrell state with Q ^ 0, the remainder domain 
is occupied by the unpaired state (A = 0). The three phases 
meet at the (Lifshitz) tricritical point. 



and for low T the FF state occupies approximately 1/4 
of the range of 8^/2 values that admit pair correlations; 
this range is reduced to 0.14 for / = 2.93. The zero- 
temperature gap is A = 64.7 MeV in weak coupling; the 
critical value of the mismatch at which the transition to 
the BCS takes place is read off to be <fyti/A(0) ~ 0.8; the 
transition from the BCS to the FF state takes place at 
<WA(0) ~ 1. 



B. Temperature-asymmetry dependence of the gap 

The detailed features observed in the phase diagram 
can be understood by examining the behavior of the gap 
function along vertical cuts of constant asymmetry and 
horizontal cuts of constant temperature in Fig. [TJ We 
shall first discuss the weak-coupling domain and then 
simply show that the behavior of the gap function in 
strong coupling is self-similar to that of weak coupling. 
Fig. [3] shows the weak-coupling gap as a function of 
temperature for a range of asymmetries. The curves 
for non-zero asymmetries can be divided into two seg- 
ments which are separated by the point where the dot- 
ted curves branch off. The high-temperature the segment 
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FIG. 3: (Color online) Dependence of the pairing gap for BCS 
and FF states on temperature for several asymmetries. The 
FF phase exists in the temperature domain between T = 
and the temperature where the dotted (red online) curves 
merge with the others. For higher temperatures the super- 
conductor is in the BCS state. The low-temperature behav- 
ior of the BCS phase is anomalous as the gap decreases with 
decreasing temperature and for asymmetries a ~ 0.07 gives 
rise to a second critical temperature. 



corresponds to the BCS state; the temperature depen- 
dence of the gap is standard, i.e., dA(T)/dT < and the 
high-temperature asymptotics is A(ev) ~ [T c (a)(T c (a) — 
T)] 1 / 2 , where T c (a) is the (upper) critical temperature. 
On the low-temperature segment there are two compet- 
ing phases (BCS and FF), which have very different tem- 
perature dependence of the gap function. The quench- 
ing of the BCS gap (dotted lines) as the temperature 
is decreased is caused by the loss of coherence among 
the quasiparticles as a consequence of vanishing smear- 
ing of the Fermi surfaces with temperature, which scales 
as [(3Ep(u, Consequently, in the low-temperature 

segment dA(T)/dT > and for large enough asymme- 
tries there exists a lower critical temperature T* [54| . 
On contrary, for the FF phase dA(T)/dT < as is 
the case in the ordinary (symmetrical) BCS theory. It 
should be mentioned that the "anomalous" behavior of 
the BCS gap below the point of bifurcation of the FF 
state leads to a number of anomalies in the thermody- 
namic quantities, such as negative superfluid density or 
negative difference between the normal and superfluid 
entropies [5J|- These anomalies are absent in the FF 
state [III. In Fig. 0] we show the gap function along cuts 



FIG. 4: (Color online) Dependence of the pairing gap for 
BCS (dashed lines) and FF (solid lines) phases on asymmetry 
parameter for several temperature labeled in the figure. 




FIG. 5: (Color online) Dependence of the gap function in 
the the FF (low temperatures) and BCS (high temperatures) 
phases on temperature for several fixed asymmetries labeled 
in figure. Left panel corresponds to the weak coupling limit 
/ = 1, the right panel to the strong coupling limit / = 2.93. 
The energy scales are given in units of the respective gaps in 
the zero temperature and symmetric limits Aoo- 



of temperature in the phase diagram. As above, there 
are two segments for each temperature: the low-a do- 
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FIG. 6: (Color online) Dependence of occupation numbers on the mode. The majority and minority components are shown for 
each asymmetry with same style curves. The numeric labels define the density asymmetry parameter a. The left, middle and 
right panels correspond to the weak f — 1, intermediate / = 2 and strong / = 2.93 coupling limits. The upper row displays 
the dependence at low temperatures, the lower row - at high temperatures (see the labeling in the figure). 



main where only the BCS phase exists and the large-cv 
domain where both BCS (dashed lines) and FF (solid 
lines) solutions exist, but the FF solution is favored. For 
small a the gap function is linear in ev; the large-a asymp- 

1/2 

totics is A(a) ~ A o (1 — a/a\) ' , where ai ~ Aqo//Z 
and Aoo is the value of the gap at vanishing temperature 
and asymmetry. The critical asymmetry ot2 at which 
the FF phase transforms into the normal phase is a de- 
creasing function of temperature, whereas that for the 
BCS phase (a± above) increases up to the temperature 
where a\ = a^} for larger temperature a,\ decreases with 
temperature. Consequently, in the dominant phase the 
critical asymmetry always decreases with temperature. 

Figure [5] compares the temperature dependence of the 
gap in the preferred phase (FF phase at low-T, BCS 
phase at high-T) in the weak-coupling (/ = 1, left panel) 
and strong-coupling limits (/ = 2.93, right panel). We 
observe that the curves in both limits are self-similar, 
therefore the qualitative arguments presented for the case 
of weak coupling should hold also in the case of strong 



coupling. The two limits differ, of course, quantitatively, 
which is manifest in the magnitudes of the asymmetries 
required to suppress the gap to a given fraction of the 
gap in the symmetric case. 



C. Occupation numbers and wave functions 

The integrands of Eqs. (|43|) and ([44]) define the occu- 
pation numbers n d / u (p) of the down and up quarks (re- 
ferred also as the majority and the minority components). 
Their dependence on momentum for a special direction 
cos 6 = is shown in Fig. [6] The weak-coupling and low- 
temperature case (left upper panel) develops a "breach" 
or "blocking region" for large asymmetries, i.e., the mi- 
nority component is entirely expelled from the blocking 
region (n u = 0), while the majority component is maxi- 
mally occupied (rid = !)• hi the small a limit occupation 
numbers are clearly fermionic (with some diffuseness due 
to the temperature), i.e., the states are filled below a cer- 



12 




P [fin"'] P [fin" 1 ] 



FIG. 7: (Color online) Dependence of occupation numbers on 
the momentum for weak coupling / = 1 and two angles 6 = 
45° (upper panel) and 9 = 90° (upper panel). The majority 
and minority components are shown for each asymmetry with 
curves of the same style. The numeric labels indicate the 
density asymmetry parameter a. 



tain mode (Fermi momentum) and are empty above. In 
the high-temperature limit (lower left corner) the breach 
gets filled in, the occupation numbers are smooth func- 
tions of momentum, consequently the high-momentum 
modes are less populated. In the intermediate-coupling, 
low-temperature, regime (/ = 2, middle panel in Fig. 
the fermionic nature of the occupation numbers is lost, 
e.g., the zero-momentum modes are not fully populated 
until large asymmetries are reached. Furthermore, it is 
not possible to identify a Fermi surface as the modes are 
populated smoothly (an exception is the population of 
the majority at near critical asymmetries). Temperature 
smears out the "knee" in the occupation numbers, as 
seen in the plot for the intermediate-coupling and high- 
temperature regime (lower middle panel). The strong- 
coupling regime / = 2.93 can be identified with the Bose- 
Einstein condensate phase of strongly coupled pairs. At 
large asymmetries the distribution of the minority com- 
ponent undergoes a topological change, by first devel- 
oping an empty strip within the distribution function, 
which is followed at larger asymmetries by a distribution 
in which modes are populated starting from a certain 
nonzero value. Thus, the Fermi sphere of the minor- 
ity component in the BCS limit evolves into a shallow 
shell structure in the strongly coupled Bose-Einstein con- 
densed limit. Such a behavior has been established both 
for nonrelativistic [4fJ and relativistic systems [44[ in the 



FIG. 8: (Color onlne) Same as in Fig. [7] but for strong cou- 
pling, / = 2.93. 



case of a homogeneous superconducting phase. 

To understand the behavior of the occupation num- 
ber in the FF phase, we need to look at the cases where 
cos 9^0 (the limit cos 6 — discussed above is spe- 
cial as it describes the asymmetric BCS state). Figure [7] 
shows the occupation numbers in the weak-coupling limit 
for two fixed angles 9 = 45° (upper panel) and 9 = 90° 
(lower panel) and for several asymmetries. For fixed an- 
gle 9 = 45° increasing the asymmetry has an effect on 
the occupation number which is opposite to the one ob- 
served for 9 — 0° case: for large asymmetries the dif- 
ference between the occupation numbers disappears, i.e., 
the superconductor behaves as if it was flavor-symmetric. 
Thus, asymmetries and angles for which pairing is pro- 
moted over the asymmetrical BCS state makes it possible 
for the FF state to have lower free energy than the BCS 
state; in effect, the nonzero momentum compensates for 
the mismatch of the Fermi spheres and restores the co- 
herence needed for pairing to occur. For fixed 9 = 90° 
the picture is reversed again: increasing the asymmetry a 
blocking region (breach) develops in the occupation num- 
bers; note however that there is no perfect symmetry in 
the occupation numbers against rotations by 90°. 

Figure [S] shows the same dependence in the case of 
strong coupling / = 2.93. Both for 9 = 45° (upper panel) 
and for 9 = 90° (lower panel) the occupation numbers 
move apart with increasing asymmetry, as is the case for 
9 = 0°. For large asymmetries and nonzero angles the 
knee in the majority distribution is less pronounced and 
the minority component does not undergo the topolog- 
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FIG. 9: (Color online) Dependence of the kernel K(p) on the momentum in the two extremes of weak / = 1 (left column) 
and strong couplings / = 2.93 (right column). The numeric labels define the density asymmetry parameter a. The upper row 
corresponds to low temperatures, the lower row - to high temperatures (see the labeling in the figure). 



ical change seen for 9 — 0°, but the general features of 
occupation numbers are the same. 

Fig. [9] displays the kernel of Eq. (I42p . which can be 
written as 



K( P ) = J2 



1 



tanh 



f£ e + (A e ) 



A e | 2 
- tanh 



f£ e "(A e ) 



.(50) 



For contact interactions this kernel is the product of 
the imaginary part of the (on-shell) retarded propagator 
and the Pauli operator (the sum of tanh terms). Physi- 
cally, it has been interpreted as the wave function of the 
Cooper pairs, since in the non-relativistic limit it obeys 
a Schrodinger-type eigenvalue equation. The prefactor 
of the Pauli operator is a smooth function of momentum 
with a maximum at the Fermi surface, where Es, e = 0. 

The modes that essentially contribute to the pairing 
correlations in different regimes can be identified from 
Fig. [9l In the low-temperature regime K(p) has two 
maxima which are separated by a depression around the 
Fermi momentum. This behavior is the consequence 
of the blocking region in the occupation number dis- 
cussed above. To see this we note that the Pauli op- 
erator (up to an irrelevant constant) can be rewritten as 
1 - fu(p) - fd(p), where f u / d (p) are the u and d quark 
distributions. We have seen that in the blocking re- 
gion fd(p) — > 1 and fuip) — * 0, as a consequence the 
function K (p) is depressed for modes within the block- 
ing region. Increasing the temperature smears out the 
structures seen in the low-temperature case, in agree- 



ment with the smearing seen in the occupation numbers. 
In the strong-coupling and low-temperature limit (left 
upper panel) the double-structure is broadened due to 
the correlations rather than the temperature. For large 
asymmetries there is a topological change in the struc- 
ture of the Fermi sphere of the minority which is seen 
to depress the low-momentum structure from the minor- 
ity component, i.e., the contribution to the pairing comes 
only from the vicinity of the Fermi sphere of the majority 
component. 

The quantity K(p)c&n be interpreted as a wave func- 
tion of the Cooper pairs. It is seen that in the weak- 
coupling limit it is well localized in momentum space, 
i.e., the Cooper pairs have structure which is broad in real 
space (large coherence length). The picture is reversed 
in the strong coupling limit, where K{p) is a broad func- 
tion of momentum, therefore it is well localized in the 
real space, i.e., it corresponds to strongly bound Bose 
molecules. The temperature suppression of the pairing 
correlation in either limit can be seen by comparing the 
upper and lower panels of Fig. 

We consider next the Fourier transform of kernel K (p) 
to understand the spatial structure of the Cooper pairs. 
The Fourier image of K{p) function has the physical in- 
terpretation of the Cooper pair "wave-function" , which 
we define as 

V-(r) = y/N [ £P [K(p, A) - K(p, 0)] e^ r 



2irr 2 



(2tt) 3 

/ dpp [K(p,A)-K(p,0)]sm(pr), 
Jo 



(51) 
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FIG. 10: (Color online) Upper panel: The wave-function of Cooper pairs for weak / = 1 (left column), intermediate / = 2 
(middle column) and strong / = 2.93 (right column) coupling. The asymmetries are a — 0.1 (solid line, black online), 0.2 
(dashed line, red online), 0.3 (short-dashed line, blue online), 0.4 (dashed-dotted line, orange online), and 0.5 (dashed-double- 
dotted line, green online). The dotted curve shows the results for a = 0.01 at intermediate coupling; in the case of weak 
coupling the a — 0.01 curve coincides with a = 0.1 one. The temperature is T — 1 MeV in the left panel, and T = 40 MeV in 
the middle and right panels. 



where N is a constant that is determined from the nor- 
malization condition 

N J d 3 r\^(r)\ 2 = 1. (52) 

In Eq. (|5ip we subtracted from the kernel its value 
in the normal state K(p, 0) to regularize the integral, 
which is otherwise divergent. Cut-off regularization of 
this strongly oscillating integral is not appropriate. The 
mean-square radius of a Cooper pair is defined via the 
second-moment integral of the density probability 

(r 2 ) = J d 3 rr 2 \Tjj(r)\ 2 . (53) 

The coherence length, i.e., the spatial extension of a 
Cooper pair, is defined as £ = \J (r 2 ). The regimes of 
strong and weak coupling can be identified by comparing 
the coherence length to the mean interparticle distance, 
which is defined as d — [3/47r(n tl +rif,)] 1 / 3 . We obtain d = 
0.6 fm, for n B = 0.55 fm~ 3 , and £(/ = 1) = 5.416942, 
£(/ = 2) = 1.9, and £(/ = 2.93) = 1.9. 



Fig. [TU1 upper panel, shows the wave-function of 
Cooper pairs as a function of radial distance for various 
coupling constants. In weak coupling, the wave func- 
tion has well-defined oscillatory form which extends over 
many periods of the interparticle distance. Such a state 
corresponds to the BCS theory where the spatial correla- 
tions are characterized by scales that are much larger 
than the interparticle distance. For intermediate and 
strong coupling the wave function is increasingly con- 
centrated at the origin with few periods of oscillations 
at most. The strong-coupling limit corresponds to pairs 
that are well localized in space within a small radius. 
This regime has clearly BEC character, where the pair- 
correlations extend over distances comparable to the in- 
terparticle distance. It seen that in weak coupling the 
wave function is almost independent of the asymmetry, 
whereas in strong coupling this dependence is substan- 
tial. The complementary lower panel of Fig. [TO] shows 
the quantity r 2 1 -0 (t) | 2 . In strong coupling the spatial 
correlation is dominated by a single peak, which indi- 
cates a tightly bound state close to the origin. The exis- 
tence of residual oscillations indicates that there is no 
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unique bound state formed at such coupling, but the 
tendency towards its formation is clearly seen. An os- 
cillatory structure appears in the intermediate-coupling 
regime indicating a transition from the BEC to the BCS 
regime. In the weak-coupling limit the short-range cor- 
relations disappear leaving behind a BCS state that is 
correlated over large distances. At low and high asymme- 
tries the strong-coupling peaks are well defined, whereas 
at intermediate asymmetries the weight of the function 
is distributed among several peaks. 



IV. CONCLUSIONS AND OUTLOOK 

We have calculated the phase diagram of chiral quark 
matter in the temperature-flavor asymmetry plane, in 
the case where there are three competing phases - the 
homogeneous BCS state, the inhomogeneous FF state 
and unpaired quark matter. The theory is formulated in 
the chiral limit of massless up and down quarks interact- 
ing with a contact interaction obtained from integrating 
out the gluonic degrees of freedom. Fluctuations beyond 
mean-field have not been included. We have studied how 
the structure of the phase diagram changes as the attrac- 
tive coupling strength is increased. At weak coupling 
and zero temperature the FF state is favored for arbi- 
trary asymmetries. At higher temperatures its domain 
of occupancy is restricted to increasingly larger asym- 
metries. The phase diagram features a tricritical Lif- 
shitz point, where the BCS, FF and normal states meet. 
Near this point our system belongs to the universality 
class of paramagnetic-ferromagnetic-helical systems, for 
which the critical exponents near the Lifshitz point are 
known. In particular, the momentum scale characteriz- 
ing the breaking of translational invariance has a critical 
exponent of 1/2 to leading order. As the coupling is in- 
creased the domain occupied by the BCS phase relative 
to the FF phase increases. There exists a critical cou- 
pling at which the FF phase disappears entirely. This 
transformation with coupling shifts the tri-critical point 
to larger asymmetries, where it disappears when the FF 
state is extinct. The strong-coupling limit of the FF state 
is particularly interesting, since at such coupling the sys- 
tem enters the BEC regime, as is evidenced by the mo- 
mentum dependence of the occupation numbers and the 
ratio of the coherence length to the interparticle distance 
being of the order of unity. This new regime can be iden- 
tified with a current-carrying Bose-Einstein condensate 
of bound chiral up and down quarks. 

The analysis of the dependence of the gap function as a 
function of temperature and asymmetry shows that this 
functional dependence carries over to the intermediate- 
and strong-coupling regimes in a self-similar manner. In 
particular, the temperature dependence of the gap in the 
BCS and FF states does not show anomalies present in 
the asymmetric BCS state, i.e., their behavior is self- 
similar to the one observed in flavor-symmetric BCS 
phase. Our analysis of the momentum dependence of 



the occupation numbers shows that in the BCS state, 
and in the FF at special angle for which the condensate 
momentum Q = 0, the minority component contains a 
blocking region (breach) around the Fermi sphere in the 
weak-coupling limit, which engulfs more low- momentum 
modes as the coupling increases, and eventually leads to a 
topological change in strong coupling, where the minor- 
ity Fermi sphere contains either two occupied strips or 
an empty sphere [4(| H] . We show that at arbitrary an- 
gles, i.e., for (effectively) non-zero Q, the blocking region 
is filled in and the momentum dependence of the occu- 
pation numbers is smooth; this anti-blocking effect is a 
direct consequence of the mechanism of restoration of the 
pair correlation in the FF state due to non-zero Q. This 
mechanism can be seen in the momentum dependence 
of the Cooper pair wave function. In the asymmetric 
BCS state the blocking of minority occupation numbers 
leads to a depletion of in the wave-function, and hence 
reduction of the gap. The anti-blocking effect in the oc- 
cupation numbers of the FF phase, inverts this tendency 
by restoring the wave- function, and therefore pair corre- 
lations. In weak coupling the momentum dependence of 
the occupation numbers is a Fermi distribution (neglect- 
ing for the sake of argument the changes near the Fermi 
surface due to pairing and flavor asymmetry). In strong 
coupling, the momentum dependence does not show a 
step- function behavior, rather one deals with a broad dis- 
tribution with modes occupied far away from the actual 
Fermi sphere. This form of occupation numbers suggest 
that we work effectively in the Bose-Einstein-condensed 
limit. This is further confirmed by our study of the 
spatial structure of the Cooper wave function. In the 
weak-coupling limit it describes a macroscopically coher- 
ent state characterized by wave function oscillations that 
extend over distances which are much larger than the 
interparticle distance. In the strong-coupling limit the 
nodal structure of the wave function disappears; instead 
one observes a well developed single structure, which is 
interpreted as a bound state of up and down quarks (and 
anti-quarks) forming a boson. This picture lends sup- 
port to the claim that in strong coupling we observe a 
new regime of the FF state which corresponds to the 
current-carrying Bose condensate. 

Minimal extensions of the present model to account for 
electric and color charge neutralities and /3-equilibrium 
among quarks will make it suitable for applications to 
compact star physics [25j, |79|, ISO, [8l| . 

The strong-coupling limit of the FF phase is not likely 
to be realized in quark matter in compact stars. Cur- 
rently accepted parameters predict quark superconduc- 
tivity well in the BCS regime. Nevertheless, our study 
of the strong-coupling limit, and identification of a new 
regime of current-carrying Bose condensate, may not be 
entirely of academic interest. Cold fermionic atoms offer 
an excellent laboratory where the properties of the phases 
of imbalanced fermions are being tested. The diversity of 
such systems and the methods (ranging, e.g., from direct 
imaging of the atomic clouds to studying the vortex lat- 
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tices) hold the promise that our findings could be tested 
in experiments [M M, HH, [H, M, M, HE HE M, H|. 
The relativistic (chiral) aspects of our study could be 
relevant for the studies of graphene - another condensed 
matter system, where electrons are described by the 
Dirac equation and, therefore, exhibit relativistic dynam- 
ics [H M m, m . 
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APPENDIX A: THE QUASIPARTICLE SPECTRUM 

In this appendix we calculate the explicit form of the quasiparticle spectrum in the FF state from Eq. (|2!l 



(W)~ l (^/)" l -Ei A "«i aAe w»= ' 



(Al) 



where the Green's functions in the chiral limit are given by 



The explicit form of Eq. (|A1|) is 



■~2 +k >* ) 7o 



£|A e ^| 2 A e (fc) = 0. 



(A2) 



(A3) 



Upon multiplying this equation from left and right by 7 one finds 



tV 1-^ + *, 



7Vf% + *, 



^|A^| 2 7 °A e (fc)7°-0 



(A4) 



We further separate the time and space components and use the projectors to positive and negative energy states 
A ± (fc) = (1 + a ■ k/\k\)/2 and their property 



|fe|[A+(fc) - A-(fc)] = a-fc, 



(A5) 



to eliminate the quantity a = 7°7. Since the dispersion relations for the particles and anti-particles separate, we 
shall keep only the particle states below to obtain 



1 



A + (k-Q) + A+\k\ + k Q - f i{ 



--A+(fe-Q)-A+|fc|+fc + M? 



|A%| 2 A+(fc) = 0. 



The averaged and mismatched chemical potentials are defined as 

l f l f 

i.e., /i{ = p, + S/j, and jjp, = Jl — Sfi. Substituting these relations we obtain 



(A6) 



(A7) 



--A+(fc • Q) + A+|fe| + k - p. - S/j, 



--A+(fe • Q) - A+|fe| + fc + ft - Sfi 



|A%| 2 A+(fc) = 0. 



(A8) 
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The terms that do not contain the positive projector can be multiplied by A + + A 
to the dispersions of the positive and negative energy state particles. We obtain 



1, i.e. they equally contribute 



k --(k-Q) 



Sfi + \k\ — fx 



ko-^{k-Q)-8ft-(\k\-fi) 



|A /! ?,| 2 = 0, 



(A9) 



or 



k — -(k ■ Q) — 5(i 



(\k\-(i) 2 -\A 



fg |2 



(A10) 



The dispersion relation is the solution of Eq. (|A10|) 



k Q = -(k-Q) + 5fi±^(\k\-fiy + \A } 



(All) 
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